Draft version February 2, 2008 

Preprint typeset using I4TgX style emulateapj v. 10/10/03 



THE SURPRISINGLY STEEP MASS PROFILE OF ABELL 1689, FROM A LENSING ANALYSIS OF SUBARU IMAGES' 

Tom BROADHURST^, MASAHIRO TAKADa\ KEIICHI UMETSU'*, XU KONg\ NOBUO ARIMOTO^ Masashi Chiba^, Toshifumi 

FUTAMASE^ 
Draft version February 2, 2008 



o 
o 

(N 
o 

CD 

Q 

00 



> 

(N 
O 



Oh 
I 

O 



X 



ABSTRACT 

Subaru observations of A1689 (z = 0.183) are used to derive an accurate, model-independent mass profile for 
the entire cluster, r ^ 2Mpc//i, by combining magnification bias and distortion measurements. The projected 
mass profile steepens quickly with increasing radius, falling away to zero at r ^ LOMpc/Zi, well short of the 
anticipated virial radius. Our profile accurately matches onto the inner profile, r ^ 200kpc//z, derived from 
deep HST/ACS images. The combined ACS and Subaru information is well fitted by an NEW profile with 
virial mass, (1.93 ±0.20) x lO'^M©, and surprisingly high concentration, Cvir = 13.7;|;] |, significantly larger 
than theoretically expected (cyir — 4), corresponding to a relatively steep overall profile. A slightly better fit is 
achieved with a steep power-law model, d\ogY,{9) / d\og9 ~ -3, with a core 6c ~ 1.'7 (r^ ~ 210kpc//!), whereas 
an isothermal profile is strongly rejected. These results are based on a reliable sample of background galaxies 
selected to be redder than the cluster E/SO sequence. By including the faint blue galaxy population a much 
smaller distortion signal is found, demonstrating that blue cluster members significantly dilute the true signal 
for r ^ 400kpc//i. This contamination is likely to affect most weak lensing results to date. 

Subject headings: cosmology: observations - gravitational lensing - galaxies: clusters: individual(Abell 1689) 



1. EvlTRODUCTION 

Numerical simulations based on the cold dark matter 
(CDM) scenario are reliable enough to make statistical pre- 
dictions for the mass profiles of clusters. The gradient of an 
'NEW profile is predicted to monotonically steepen with in- 
creasing radius (Navarro, Frenk & White 1997), with loga- 
rithmic slopes shallower than an isothermal profile interior 
to the characteristic radius r < r,, but steeper at larger ra- 
dius, approaching to at r ^ ryh-. This curvature is par- 
ticularly pronounced for massive clusters, where the halo is 
expected to have a relatively low concentration, providing 
a clear prediction. Weak lensing work has yet to make a 
definitive statement regarding the validity of the NEW pro- 
file, with the analysis claiming only broad consistency with 
both the singular isothermal case and the NEW model (Clowe 
& Schneider 2001; Bardeau et al. 2004). Recently, Kneib et 
al. (2003) made a more comprehensive lensing study using 
the panoramic sparse-sampled HSTAVFPC2 data of C10024, 
showing that at large radius, the mass profile is steeper than 
the isothermal case, consistent with an NEW profile for r> r,. 
A firmer constraint is possible with magnification informa- 
tion, breaking the mass-sheet degeneracy inherent to dis- 
tortion measurements (Broadhurst, Taylor & Peacock 1995, 
hereafter BTP). Detections of a number-count depletion have 
been claimed for several clusters (Fort et al. 1997; Taylor et 
al. 1998; Athreya et al. 2002; Dye et al. 2002). 

Recently, the central region of A 1689 (z = 0.183) has been 
imaged in detail with the Advanced Camera for Surveys 
(ACS), revealing 106 multiply lensed images of 30 back- 
ground galaxies (Broadhurst et al. 2004, B04). Many radially 
directed images define a radial critical curve, inside which 
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small counter-images are identified, so the mass profile can 
be traced in detail to the center of mass. An NEW profile fits 
well over the restricted ACS field, r ^ 200kpc//!, but with a 
somewhat larger concentration, Cyh = 8.2^j g, than expected 
for massive clusters, Cyir ~ 4 (e.g.. Bullock et al. 2001). To 
examine the shape of the full profile, we turn to the wide-field 
prime focus camera, Suprime-Cam, of the 8.2m Subaru tele- 
scope. Suprime-Cam provides an unparalleled combination 
of area (34' x 27') and depth (Miyazaki et al. 2002). We mea- 
sure the lensing distortion and magnification of background 
red galaxies and combine them to derive a model-independent 
mass profile out to r ~ 2Mpc, allowing a definitive compari- 
son with model profiles. Throughout this Letter, the concor- 
dance ACDM cosmology is adopted (fimo = 0.3, JIao = 0.7, 
h = 0.7). Note that one arcminute corresponds to the physical 
scale 129kpc//i for this cluster 

2. DATA REDUCTION AND SAMPLE SELECTION 

Suprime-Cam imaging data of A 1689 in V (1,920s) and 
SDSS /' (2,640s) were retrieved from the Subaru archive, 
SMOKA. Reduction softwai-e developed by Yagi et al (2002) 
was used for flat-fielding, instrumental distortion correction, 
differential refraction, PSE matching, sky subtraction and 
stacking. The resulting FWHM is 0".82 in V and 0".88 in 
/' with 0".202 pix"', covering a field of 30' x 25'. 

Photometry is based on a combined V + Z' image using SEx- 
tractor (Bertin & Arnaut 1996). The limiting magnitudes are 
V = 26.5 and /' = 25.9 for a 3cr detection within a 2"aperture 
(we use the AB magnitude system). For the number counts to 
measure magnification, we define a sample of 8,907 galaxies 
(12.0 arcmin"^) withV-/' > 1.0. For distortion measurement, 
we define a sample of 5,729 galaxies (7.59 arcmin"^) with 
colors 0.22 mag redder than the color-magnitude sequence 
of cluster E/SO galaxies, (y-/') + 0.0209/'- 1.255 > 0.22. 
The smaller sample is due to the fact that distortion analy- 
sis requires galaxies used are well resolved to make reliable 
shape measurement. We adopt a limit of /' < 25.5 to avoid 
incompleteness effect. Our red galaxies are very reasonably 
expected to lie behind the cluster, made redder by larger k- 
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corrections. The counts of all galaxies including the dominant 
faint blue population is much larger, ^ 40 arcmin"^. How- 
ever, as we show below, the distortion signal of the full sam- 
ple is significantly diluted within r ^ 400kpc//!, compared to 
the red sample, indicating that blue cluster members contam- 
inate the sample. The mean redshift of the red galaxies is 
estimated to be ~ 1 ±0.1, based on deep photo-z estima- 
tion for deep field data (Benitez et al. 2002). In the following 
we will assume (zj) = 1 for the mean redshift, but note that the 
low redshift of A 1689 means that for lensing work, a precise 
knowledge of this redshift is not critical. 

3. LENSING DISTORTIONS 

We use the IMCAT package developed by N. Kaiser for our 
distortion analysis, following the well tested formalism out- 
lined in Kaiser, Squires & Broadhurst (1995), with modifica- 
tions described by Erben et al.(2001). 

Briefly, to estimate the lensing distortion from the ob- 
served ellipticities of galaxy images, we first correct for PSF 
anisotropy using the available stars. The mean stellar ellip- 
ticity before correction is (e"i*,e"2*) ~ (-0.013,-0.018) over 
the survey field, and is reduced by the correction to only 
= (0.47 ± 1.32) x lO'^ and e*""' = (0.54 ±0.94) x lO'^^. 
We also correct for the isotropic smearing effect caused by 
seeing as well as by the window function used in the shape es- 
timate. Full details will be presented in Umetsu et al. (2004, 
in preparation) including the 2D maps. 

Figure [2 shows the radial profiles of the tangential distor- 
tion, and the 45-degree rotated (x) component, gx, for 
the different samples of galaxies, where the cluster center is 
well determined by the locations of brightest cluster galaxies 
and/or giant arcs. Note that lensing induces only the tangen- 
tial shear in the weak lensing regime. The error bars are 68% 
confidence intervals computed from the bootstrap resampled 
data. It is clear that the observed signal is significant to the 
limit of our data, 20' or 2Mpc / h, and that the x -component 
of the red galaxy sample is consistent with a null signal at all 
radii, indicating the reliability of our distortion analysis. 

One can see that the NFW prediction for source redshift 
z, = 1 (cvh- = 8.2 and Myir = 2.6 x lO'^/iM©), which fits best 
the ACS strong lensing data restricted to the central region 
^ 2' (B04), matches onto the Subaru distortions at ^ 3' for 
the background red galaxies, but increasingly overestimates 
the distortions at larger radii ^ 3'. A steeper NFW profile 
with Cvir ~ 14 (solid curve) better reproduces the distortion 
profile at all radii, as does a cored power-law profile (dotted 
curve), (see ^for the details). 

A careful background selection is critical. If we select all 
galaxies irrespective of color with magnitude limit /' < 25.5, 
and exclude only the cluster sequence galaxies, then we find 
that the distortion signal falls below the red background sam- 
ple by a factor of 2-5 at r ;S 400kpc//i. While the slight dilu- 
tion within a factor of 2 is still apparent at large radii ^ 5', 
blue cluster members must be significantly contaminating the 
full sample at the small radii, reducing the distortion signal. 
This full sample selection is very similar to Bardeau et al. 
(2004) and Clowe & Schneider (2001), and we find close 
agreement. Hence it is apparent that these analyses under- 
estimate the true distortion signal and explains why they un- 
derpredicted the Einstein radius, 20", compared with the 
much larger observed radius of ^ 45" (for Zs = 1), a puzzle 
noted by Clowe & Schneider (2001). 

It is also instructive to examine the tangential distortions of 
galaxies whose photometry places them in the limits bound- 
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Fig. 1 . — Radial profiles of the tangential shear (upper panel) and the 45° 
rotated (x) component (lower panel). The square symbols correspond to the 
red background sample. The cross symbols con'espond to the full sample of 
faint red+blue galaxies, excluding objects falling on the narrow cluster E/SO 
sequence. The distortion signal of the full sample falls well below that of the 
background red galaxies for r ^ 400kpc / h, indicating significant contami- 
nation by blue cluster members. The triangle symbols denote galaxies with 
photometry coincident with the color-magnitude sequence of cluster E/SO 
galaxies (the open symbol indicates the negative signal). The solid and dot- 
ted curves show an NFW profile with a large concentration, fyii ~ 14, and 
a power-law profile with slope « = 3.16 and core 9c = 1.65 arcmin (see jsl . 
respectively, both matching well the overall profile of the red galaxy sample. 
The best-fitting model (dashed curve) based on the inner ACS measurements 
(B04), with c'vii ~ 8, clearly overestimates the outer profile. 



ing the cluster sequence. A null signal is found for ^ 5' as 
expected, where the cluster member is prominent above the 
background. But a positive distortion, with weaker signal than 
the red galaxy sample, is found for ^ 5', indicating that back- 
ground, lensed galaxies contribute significantly. 

4. MAGNIFICATION BIAS 

Lensing magnification, /i(0), influences the observed sur- 
face density of background galaxies, expanding the area of 
sky, and enhancing the flux of galaxies (BTP). The num- 
ber counts for a given magnitude cutoff m^ut, approximated 
as a power law with slope s = d\ogN(,{m) / dm, are modified 
as A^(< mcut) = Nq{< OTcut)/i^'^'~', where A^o is the unlensed 
counts. Thanks to the large Subaru FOV, the normalization 
and slope of A^o for our red galaxy sample are reliably esti- 
mated as «() = 12.6 ± 0.23 arcmin"^ and s = 0.22 ± 0.03 from 
the outer region > 10'. The slope is less than the lensing in- 
variant slope, s = 0.4, so a net deficit of background galaxies 
is expected. 

We conservatively account for the masking of observed sky 
by excluding a generous area irab around each cluster se- 
quence galaxy, where a and b are defined as 3 times the ma- 
jor (A_IMAGE) and minor axes (B_IMAGE) computed from 



Broadhurst, Takada, Umetsu et al. 



3 



200 



r [kpc//i] 
500 1000 



2000 



3 




ACS+Subaru NFW profile 
ACS+Subaru CPL profile 

■ ACS NFVr profile 

■ red galaxies (masking corrected) 
A red galaxies (no masking) 



5 10 
6 (arcmin) 



20 



Fig. 2. — Number-count profile of the red background galaxies (square 
symbols). The counts are normalized by the unlensed counts estimated from 
the cluster outer region (see text for details). The triangle symbols show 
the counts without the mask correction due to cluster members. The model 
curves are shown for comparison as in FigurelTI 



SExtractor, corresponding to the isophotal detection limit 
(Berlin & Arnouts 1996). The number density in each radial 
annulus is then calculated by excluding this area and simply 
renormalising. Note that, if we adopt the masking factor of 2 
or 4 instead of 3, the results shown below are little changed. 

Figure 121 shows that the red galaxy counts are clearly de- 
pleted, with a clear trend towards higher magnification in the 
center The masking area is negligible at large radius and rises 
to ^ 20% of the sky close to the center, r ^ 3'. Comparison 
with models shows that the data are broadly in agreement with 
the profiles as described in |5] We have ignored the intrinsic 
clustering of background galaxies, which seems a good ap- 
proximation, though some variance is apparent in 2D maps 
(Umetsu et al. 2004) and may explain a discrepant point at 
~ 5' in the radial profile. 

5. MODEL-INDEPENDENT MASS PROFILE 

The relation between distortion and convergence is non- 
local, and masses derived from distortion data alone suf- 
fers from a mass sheet degeneracy. However, by combin- 
ing the distortion and magnification measurements the con- 
vergence can be obtained unambiguously. Here we derive a 
model-independent, discrete convergence profile in 10 log- 



arithmically spaced bins for I' < 9 < 18' 



n{9i) for 



/ = 1 , 2, . . . , 10, representing 10 free parameters which are con- 
strained with 20 data points in Figures [T]and|2l 

To perform the mass reconstruction, we need to express the 
lensing observables in terms of the binned convergence pro- 
file. The tangential shear amplitude in the /-th radial bin can 
be expressed as 7+/ = k, - k,- (Fahlman et al. 1994), where 
is the average convergence interior to radius 6*, and expressed 
as Ki = K(9i) = 2/9fJ2'j^^ 9jnjh with the bin width h = A\n9. 
Once the shear profile is given, it is straightforward to com- 
pute the binned tangential distortion and magnification via 
the relations g+i = 7+,/(1-k,) (at radii of our interest) and 
1(1 - K,)^ - 7^,1"'. However, for the discrete model, to 



the well constrained ACS derived profile. We have checked 
that, if we instead adopt the model independent mass interior 
to the Einstein radius (« 45"), the results are very similar. 

The best-fitting model k, is then derived by properly 
weighting signal-to-noise (S/N) ratios of the distortion and 
magnification data (e.g., see Schneider, King & Erben 2000). 
Since covariance between the distortion and magnification 
signals can be ignored, the is simply expressed as 

= XI + Xl with xl = E,=i [^r^^'-^Sf /4 and xl = 

E'fi [Nf°'^''^-Nf']^ /N°^'. Note that A^,- denotes the num- 
ber counts in the /-th radial annulus and the model counts 
^modei given in terms of the unlensed surface density no and 
the model magnification as N^™'^^^ = /iP'"'«()Jl,-, where il,- is 
the effective observed area of the /-th annulus excluding the 
masking area by the cluster members (see 

The square symbols in Figure |3] show the best-fitting mass 
profile obtained from the Subaru data, where the minimized 
X^ value is Xmin/'^°f- ~ 27.9/10 for 10 degrees of freedom. 
The bad Xmin mainly due to the two data points at ^ 3' 
and 5' in Figure|2l which apparently deviate from a continu- 
ous curve and are probably affected by the intrinsic clustering 
of background galaxies. Excluding the two data points leads 
to an acceptable fit x^j^/dof. = 12.4/8, meaning the best-fit 
mass profile well reproduces the distortion and magnification 
data simultaneously''. The outer profile at 9 ^9' is consistent 
with a null signal to within la (the points at ~ 9' and 16' are 
1(7 upper limits), implying a more rapid decay in mass pro- 
file than expected from the standard NFW profile, falling well 
short of the anticipated virial radius. We have also checked 
that this result is not sensitive to varying the count-slope and 
the unlensed background galaxy density within the la mea- 
surement errors. 

Force fitting an NFW profile to the Subaru mass profile 
yields Mvii- = 1.69!;jjg x IO'^Mq (r^i, = 1.95 ± 0.1 IMpc//;) 
for the virial mass (radius), where Xmin/'^'^^- - 5.36/8. Note 
that the NFW profile is specified by the two parameters, Mvi, 
and halo concentration Cvh (Bartelmann 1996), and we have 
adopted the flat prior c^h < 30 because the NFW profiles with 
Cvir ^ 20 can not be distinguished by the Subaru data alone 
due to lack of information on the inner density profile. The 
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properly compute the shear in the first bin, we need to specify 
the mass interior to radius 1', which is readily obtained from 



A pure power-law profile given by k oc 6* " also provides a 
good fit with n = 2.08^^1° (xiin = 4.07/8), corresponding to a 
3D profile, p oc r""* "^. A singular isothermal sphere with n = 1 
is excluded at more than 4a level 

Next, we consider the combined Subaru and ACS profile, 
where the ACS profile with error bars is taken from Figure 22 
in B04 and the amplitude is scaled to Zs = I from the = 3 
result. By combining the strong and weak lensing analyses, 
we can trace the mass distribution over a large range in am- 
plitude K ~ [10"^l] and in radius r = [10"^,2]Mpc//i. In 
this case, we have 22 independent data points in total exclud- 
ing the ACS data points at radii overlapping with the Subaru 
data. The best-fitting NFW profile is given by Cvir = 13.7;!;j j 
and Mvh- = (1 .93 ± 0.20) x IO'^Mq (r^i, = 2.04 ± 0.07Mpc//z), 
and has an acceptable fit Xmin/'^'^^- - 13.3/20. The inset plot 
shows how the constraints are degenerate in the (cyir, Mvii) 
plane. Hence, this result gives the lensing based confirmation 
of the NFW profile over the radii we have considered, how- 

* If we use the distortion data for the red+blue galaxies (the cross symbols 
in FigureQ, the minimum gets worse, Xmin/''°^- ~ 14-7/10, and the mass 
profile is significantly underestimated. 
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Fig. 3. — Reconstructed mass profile. The triangle and square symbols with error bars show the results from the ACS strong lensing analysis (B04) and 
the Subaru weak lensing analysis (this work), respectively. The dashed and solid curves show the best-fitting NFW profiles for the ACS data alone and for 
the combined ACS+Subaru profile, respectively. The best-fitting NFW profile for the ACS-HSubaru profile has a high concentration, Cyii = 13.7, and somewhat 
overestimates the inner slope and is a bit shallow at large radius. A cored power-law provides a better fit to the full projected mass profile. The inset plot shows 
the 68%, 95% and 99.7%. confidence levels in the (c^,„ M^„) plane for the ACS-HSubaru NFW fitting (Ax^ = 2.3, 6.17 and 11.8). 



ever, the mass distribution appears to be much more concen- 
trated toward the center than the CDM simulations predict for 
a halo of the above mass, Cvh = 4.0 (Bullock et al. 2001). A 
generalized NFW profile given by p oc r"'^(l + r/rj)"'^ (e.g., 
Moore et al. 1998) is disfavored (x^in/dof. = 28.3/20; see 
Umetsu et al. 2004 for more details), being too steep in the 
center A cored power-law profile, n cx {0 + 6c)~", gives a bet- 
ter fit: Xmin = 4.49/19 with a steep slope, n = 3.16:l:o72, and a 
core of 9, = 1.65!*];^^ arcmin (r, = 214!^^ kpc/h). A softened 
isothermal profile is strongly rejected (IQa level!). 

6. DISCUSSION AND CONCLUSIONS 

We have obtained a secure, model-independent mass pro- 
file of A 1689 over I' < 9 <20' (lOOkpc//; <r< 2Mpc//!) 
by combining the distortion and magnification-bias measure- 
ments from high-quality Subaru imaging. We have seen that 
to reliably measure distortions it is critically important to se- 
curely select background galaxies in order to avoid dilution 



of the distortion signal by blue cluster members and fore- 
grounds (see Figure [0. Thus we have resolved the discrep- 
ancy between the small Einstein radius (~ 20") inferred from 
the previous work based on largely monochromatic measure- 
ments and the observed radius 45"). The mass profile of 
A 1689 obtained from the Subaru and ACS data covers 2 or- 
ders of magnitude in radius, [10~^,2]Mpc//!, and showsacon- 
tinuously steepening profile with increasing radius, very sim- 
ilar to an NFW profile but with a much higher concentration 
than expected. The best-fitting NFW profile has Cvir — 14, sig- 
nificantly larger than expected Cyir — 4, corresponding to the 
profile expected for a much lower mass halo of ^ 1O"~'^M0 
(Bullock et al. 2001), 3-4 orders of magnitude less than the 
mass of A1689. A higher concentration may imply a higher 
than expected redshift for cluster formation, corresponding to 
a greater mean cosmological mass density, so that collapsed 
objects of a given mass have a higher internal density. A 
higher redshift of cluster formation may also help account for 
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the lack of evolution observed in the properties of galaxy clus- 
ters (Mathis, Diego & Silk 2004). Although A1689 is a very 
round shaped cluster with evidence of only modest substruc- 
ture (B04; also see Andersson & Madejski 2004 for related 
discussion based on the X-ray data), projection of structure 
along the line of sight may potentially boost the surface den- 
sity increasing the derived concentration. Substructure may 
be examined further with redshift measurements for this clus- 
ter. Simulations may also provide a good guide to the level of 
bias one may expect due to projection in general (e.g., Clowe, 
Lucia & King 2004). 

Careful lensing work on other clusters may also point to 
higher concentrations than expected. Gavazzi et al. (2003) 
quote c ~ 12 for the cluster MS2137-23, but with a low sig- 
nificance. The situation for C10024 might be complicated 
by prominent substructure, with a high concentration for the 
main clump, c ~ 22, and a low concentration, c ~ 4, for the 
secondary clump (Kneib et al. 2003). Subaru imaging in more 



passbands will provide many reliable photometric redshifts 
for refining the background selection, helping to improve the 
distortion signal and allowing the lens magnification to be 
separated from the background clustering. In addition, the 
redshift dependence of the distortion and magnification sig- 
nals may be extracted to constrain cosmological parameters. 
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